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and new physics beyond the Standard Model. Expressions for these decays are given 
in the Appendices. 
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I. INTRODUCTION 



One of the main goals of future colliders such as the CERN lhc and the future Interna- 
tional Linear Collider (ilc) is to study the Higgs sector of the Standard Model (sm) 
Moreover, the scalar sector of the SM can be enlarged with a simple extension such as 
the Minimal Supersymmetric Standard Model (mssm) and the Two Higgs Doublet Model 

nnn 

(2hdm)|4|, |a |f|, which are studied in this paper. Both in the 2hdm and MSSM the elec- 
troweak symmetry breaking is generated by two Higgs doublets fields $1,2, the imaginary 
parts of which combine to produce a CP-odd Higgs boson A and a Goldstone mode G , 
both of which must have the same quantum numbers. Consequently A and G° have the 
quantum numbers H . Owing to the J PC (A°) = + ~ the vertices A°WW and A°ZZ are 
forbidden at tree level. 

Another way to explain the absence of the A°WW and A°ZZ vertices is the following: 
these interactions come from the kinematic term after replacing one of the $'s 

by its vev, but in a CP-conserving model the VEV is real while A comes from the imaginary 
part of $°, therefore A cannot couple to any massive vector boson. Consequently these 
vertices can appear first at the one-loop level. In this study, our concern is the CP-odd A 
decay into a pair of gauge bosons, and we will therefore review the production mechanisms 
for A . 

At future e + e~ machines, the CP-even Higgs h° is mainly produced via the associated 
production with a Z boson. At high energies, WW fusion to h° becomes significant. Due 
to its CP nature, the CP-odd A possesses no tree-level coupling A°ZZ and A WW. Con- 
sequently, both of the above processes, Higgsstrahlung and WW fusion, are unavailable to 
the A at e + e~ machines. The dominant contribution is therefore from higher order dia- 



ll| and e+e~ -> A 



vv 



13j. The 



grams such as e + e — > 7 A 7, 3,3, e+e" -> ZA° l 
production rate of those loop mediated process turns out to be quite small. Therefore, the 
only realistic tree level mechanisms for CP-odd Hi ggs b osons at e + e~ colliders [1J, LUj [ly] 
are via e + e~ — > h°A° and e + e" — * bbA°, ttA° \v\. \lHAl^. At 77 collider (muon collider), A 
can be copiously produced as a resonance 77 — > A (u + u~ — >• A ) as well as in association 
with a Z boson 77 — > ZA° Q (jj+fji- ^ZA°^^i). 

The situation is different at hadron colliders. The main production mechanism is through 
gluon-gluon fusion gg — ► A or gg — > A°QQ (where Q = t. b)[2^. Those processes can 
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provide an ample sample event at the lhc. Other mechanisms like associate production 
of A with a Z boson is also possible at the lhc 24]. From the experimental side, current 
mass bounds from LEP-II for the neutral pseudoscalar A of the MSSM is m^o > 85 GeV 
in the m™ ax scenario [3, 2f|. Within 2HDM, the OPAL collaboration has used the pair- 
production e + e~ — > h°A° cross-section assuming 100% decays into hadrons independent of 
hadron flavor. It is found that the region 1 < m h < 55 GeV and 3 < m^o < 63 GeV is 
excluded at 95% CL independent of the choice of the 2HDM-II parameters |27||. 

The aim of this paper is the study of the decay of the CP-odd A into a pair of gauge 
boson both in 2hdm and MSSM. Although these decays are rare processes, loop and/or 
threshold effects could potentially give a substantial effect. Moreover, once worked out, 
any experimental deviation from the results within such a model should bring some fruitful 
information on the new physics and allow one to distinguish between models. We would like 
to mention also that these channels have a very clear signature and might emerge easily at 
future colliders. For instance, if A — > ZZ is enhanced enough, this decay may lead to the 
so-called gold-platted event signature A — > ZZ — > . 

These decays have been already studied in literature both in the mssmJ^J, S| and 
2hdm. The MSSM contribution to A — > ZZ, W + W~ has been first studied in |23j, however 
only heavy SM fermions were included. In chargino and neutralino contributions to 
A — > ZZ, W + W~ were included together with the SM fermions. However in Ref. only 
A — > ZZ was considered and the results was presented in terms of decay width only, no 
branching ratio was given. At tree level, CP-odd pseudoscalar can decay, when its kine- 
matically possible, to a pair of standard model fermions //, a pair of charginos, a pair of 
neutralinos, a pair of scalar quarks (squarks), as well as to Z°h°, Z°H° and W ± H^ . On the 
other hand, decays modes like A — >■ gg, A 77, A -> jZ, A ZZ and A W + W~ 
are mediated at one-loop level. The decays A — > gg, A — > 77, A — » 7Z has been evaluated 
by many groups and are well implemented in HDECAY program 3, [3] • In this section, we 
would like to focus on the evaluation of A — > ZZ and A —>■ W + W~ both in 2hdm and 

MSSM. 

We will update those analysis in the MSSM by including both SM fermions, charginos and 
neutralinos contributions and also by presenting both the decay width and branching ratios. 
We will also study those decays A — > ZZ, W + W~ in both type I and type II 2hdm. 

In this paper, we give a complete calculation of the CP-odd Higgs boson A decay into 
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A — > ZZ, A — > W + W~ in the general CP-conserving 2hdm and the full MSSM. In Sec- 
tion m we present the relevant Feynman diagrams and formulas involved in our calculation. 
In Section we present our numerical results and discussions. Section HXl contains our 
conclusions. 

II. CP-ODD DECAY: A -> VV 

In the MSSM or 2hdm with CP conservation, the A°ZZ and A°WW couplings are forbid- 
den at tree-level. They can be generated at the one-loop level via the Feynman diagrams in 
Fig. and Fig. EJ In Fig. topologies like 04 and 05 exactly vanishes, since Z behaves like 
CP-odd there is no transition Z-(j), <fi = h°, H°. For the same reason, A is CP-odd, there is 
no transition A°-<f) and then topologies like 06 vanishes. 

In Fig. El topologies like 08 and 09 exactly vanishes for an on-shell W boson since V4 /± -if =F 
transition is proportional to the W momentum. In topologies like 011 the transition A°-V 
is proportional to A momentum, once contracted with VW + W~ coupling it vanishes for 
an on-shell W boson. Like in Fig. 0there is no transition A°-(j) and then topologies like 010 
vanishes. 

In both cases, we did not include gauge bosons and scalars (scalar bosons and scalar 
fermions) in the loops. As it is argued in P, I3| , the sum of each of these contributions must 
cancel. The reason is that P and C are separately conserved in the bosonic sector before 
the introduction of fermions. Since the A and the Z boson are C-odd and the CP-odd 
zero angular momentum state of W + W~ must be c-even. Therefore, the coupling A°VV is 
forbidden to all order in the bosonic sector. When fermions are introduced, C, P and CP are 
no longer conserved and as a consequence A°VV can be induced at higher loop levels. 

The above statement applies to the contribution of scalar fermions as well. In this case, 
this can be easily be seen because CP-odd A couplings to scalar fermions satisfy the following 
relation A°fif* = —A°f?fj. Consequently, the total contribution of scalar fermions to A°ZZ 
or A°W + W~ exactly cancel when CP is conserved. 

Given the fact that A is CP-odd, the effective Lagrangian for A°VV must have the form: 



£>a°vv = 9a°vv^ p VfMuVps (1) 
The effective coupling gA°vv has a dimension -1, it is expected to be of the form: Qa°vv — 
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1/mw F(Ms), where F is a dimensionless function of Ms which behaves like log (Ms) and 
Ms is the masses of internal particles in the loops. Therefore, for A — > V1V2 we expect only 
a logarithmic dependence on the internal masses. The general one-loop amplitude takes the 
following form 

M(A° -> V X V 2 ) = Jt 2 ° e, vpS e^p p lP s 2 A Vl v 2 , (2) 



where Nq = 3 for quarks and 1 for leptons and charginos neutralinos in the internal loop. 
Analytical expression for Ay x y 2 both for SM fermions and charginos neutralinos loops is given 
in Appendix C. The partial decay is then computed from the above amplitude and is given 
by: 

6/\r2 \3/2 

T(A° -> v 1 v,) = s 12 J 5 C , 3 \A VlV2 \\ (3) 



where Su = 1/2 in case of identical particles in the final state and A = (m^ ~ m v 1 ~ m 
4m?/, m 2 



Both in 2hdm in MSSM, the total decay widths T^T M and T^ SM are computed as follows: 



r 2HDM = J2 T{A° -> ff) + J2 r ( A ° ^ VlV 2) + r ( A ° ^ Z$ °) + r ( A ° ^ W^ff*), ( 4 ) 
/ vs 

where the summation in Y2v ^(^4° — ► V1V2) stand for Vi — 7, g, Z and W, $° = /i° or H°. 
QCD corrections to A — * // and A — > {1/*$} decays are not included in the widths. In 
the case of MSSM, r^o SM is obtained from r^o DM by adding the decay of A® to SUSY particles: 
A - xU°j, A - xtxj and A - /J/. 

There are two types of the 2hdm, depending on which Higgs field is responsible for 
the masses of quarks and leptons. Consequently, the Yukawa couplings of the quarks and 
leptons to the Higgs bosons are different. For the A we find C = —igAOffflsf'A where the 
couplings constants are given in table [I] 

We have evaluated the one-loop induced process A — > V1V2 in the 't Hooft-Feynman 
gauge, and using dimensional regular iz at ion (3?], Q]- The types of Feynman diagrams are 



depicted in Figure Q and Figure 121 All the Fey nman diagrams have been generated and 
computed using FeynArts and FeynCalc packages. We have also used the 



e (35], l^l i: 
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(1.1) (1.2) (1.3) (1.4) (1.5) 




(1.6) 

FIG. 1: Generic contribution to A — > ZZ in 2hdm and MSSM, F denotes any fermion particles, S 
denotes any scalar particles, P denotes anay particles in the model that can fit into the diagrams 
and (f> is one of the CP-even scalars hP or H°. 



Couplings 


2HDM-I 


2HDM-II 


9A°tt 


^cot/3 


2* cot/? 


9A°bb 


^cot/3 


^tan/3 


9A0 T f 


^cot/3 


^tan/3 



TABLE I: Couplings of CP-odd Higgs boson A with fermions both in 2HDM-I and 2HDM-II. 
The MSSM is required to have 2HDM-II couplings. 

III. NUMERICS AND DISCUSSIONS 



In our numerical evaluations, we use the following experimental input quantities |3 
a^ 1 = 129, m z , m w , m t , m b = 91.1875, 80.45, 174.3, 4.7 GeV. In the MSSM, we specify 
the free parameters that will be used as fellow: i) The MSSM Higgs sector is parameterized 
by the CP-odd mass m^o and tan (3, taking into account radiative corrections from 0, Q|' 
and we assume tan (3 > 3. ii) The chargino-neutralino sector can be parameterized by the 
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P ^ 

(2.8) (2.9) 



P W 

*\cf A w 



A 

W ^ P 




(2.H) 



FIG. 2: Generic contribution to ^4° — > W + W~ in 2hdm and mssm, F denotes any fermion particles, 
S denotes any scalar particles, P denotes anay particles in the model that can fit into the diagrams 
and <b is one of the CP-even scalars h° or H°. 



gaugino-mass terms M 1; M 2 , and the Higgsino-mass term \x. For simplification M\ ~ M 2 /2 
is assumed, iii) Sfermions are characterized by a common soft-breaking sfermion mass 
Msusy = Ml = Mr, [i parameter and soft trilinear couplings for third generation scalar 
fermions A t fi, T - 

When varying the MSSM parameters, we take into account also the following constraint 
The extra contributions to the 5p parameter from the Higgs scalars should not exceed the 
current limits from precision measurements 3]: \5p\ ^ 0.003. For illustration we show first 
in Fig. El and Fig. 0] different partials decays widths of A both in 2hdm type II and type I 
as a function of m^o and for small tan j3 = 0.42. 

In Fig. El (left) we illustrate branching ratios of A as a function of m^o for small tan f3 = 
0.42 scenario in the framework of 2hdm type II. As it can be seen from the plot, for m^o < 
200 GeV, bb mode is the dominant decay mode, but when approaching ti threshold the 
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FIG. 3: Decays widths (left) and branching ratios (right) of CP-odd ^4° as a function of m^o in 2hdm 
Type II , for the parameters m h o = 130 GeV, m H o = 345 GeV, m H ± = 340 GeV, tan j3 = 0.42. 
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FIG. 4: Decays widths (left) and branching ratios (right) of CP-odd A as a function of m^o in 
2hdm model type I, with parameters same as in FigEl 

A° 9g .node is enhanced and becomes the dominant decay .node for m « £ [200,350] 
GeV. In the case when m^o « 2m t , one can say that the CP-odd is almost fermiophobic |4|. 
After crossing ti threshold, A — > ti becomes the dominant decay mode. In this scenario the 
decays of our concerns A — > W + W~ and A — > ZZ can be respectively of the order 10 _1 
and 8 x 10~ 3 near ti threshold region and decreases away from this region. 

In type I 2hdm, the situation is slightly different. For small tan/5 ps 0.42, the coupling 
of A to a pair of bottom quark, because it is proportional to cot (3, is enhanced. In this 
case, before ti threshold m A o < 2m t , bb is indisputably the dominant decay mode followed 
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FIG. 5: Branching ratios of CP-odd ^4° into ZZ, W + W~ and gg as a function of tan/3 in type II 
2hdm model (left), type I 2hdm model (right), for the parameters m h o = 130 GeV, m H o = 345 
GeV, m H ± = 340 GeV, and for various values of m^o. 

by A — > gg mode which is in the range ~ lO^-lO" 1 . Again, after crossing ti threshold, 
A — > ti becomes the dominant decay mode. As one can see, near ti threshold the branching 
ratio of A — > W + W~ and ^4° — > ZZ are of the order 10~ 2 and 10~ 3 respectively. 

In Fig. (left) we show branching ratios of ^4° as a function of tan j3 for various choice 
of rriAo both in 2hdm type II (left) and 2hdm type I (right). It is clear that in 2hdm type 
II and for small tan/3 < 1, A — > gg is the dominant decay mode for m^o < 2m t . One 
could say that in this region the CP-odd is almost fermiophobic. One can see that all the 
decays A — > VV, V = W, Z are enhanced for the small tan/3 limit. This is mainly due to 
top Yukawa coupling which is proportional to cot (3. For large tan/? all the branching ratio 
A VV decreases. 

We note finally that in 2hdm type I, the A — > gg decay mode is independent of tan f3 for 
m>Ao = 200 GeV. This is because for this value of m^o = 200 GeV, only fermionic decay are 
open and then tan f3 Dependance drop in the ratio. 
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FIG. 6: Charginos and neutralinos contributions to the decays widths and branching ratios of 
A — ► ZZ, and A — ► W + W~ as function of // parameter in MSSM and Type-II 2hdm models, for 
the parameters M Susy = 500 GeV, M 2 = 140 GeV, m A o = 350 GeV, tan/3 = 5, A t = 1 TeV. 



In the framework of MSSM, we show first in FigJ^Jthe partials decays widths FigJS](left) 
and branching ratios FigJBj (right) of A into W + W~ and ZZ as a function of \x and compare 
it to its type II 2hdm. As it can be seen from this plot, that chargino-neutralino contribution 
can either enhance the width (for Positif fi > 0) or suppress (for negative //). It is clear 
from the right plot that the branching ratio of A — > W + W~ (resp A — > ZZ) can be of the 
order 10~ 3 (resp 10" 4 ). 

In Fig. [7| and Fig. [Sj we fix \i to 1 TeV and Mi = 170 GeV and plot the decay width 
and branching ratio of A as a function of m^o for low tan/5 = 2.7 (Fig. [TJ) and tan/3 = 20 
(Fig. IBJ). For low tan (3 = 2.7, before the opening of tt mode for m^o w 350 GeV, the 
dominant decay mode for ^4° is 66 mode followed by Zh° mode for m^o > 200 GeV. Once 
tt mode is open it is fully dominating. The branching ratio of A — > WW (resp A ZZ) 
can reach a value of the order 10~ 3 (resp 10 -4 ) for CP-odd mass close to 2 m t . However for 
large tan/3 = 20, the situation is different, bb mode is fully dominating for all m^o range. 
A -> gg, Zh° are at the level of 10" 4 while A -> W + W~ and A -> are at the level of 
10 -5 and 10~ 6 branching ratio respectively. 
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FIG. 7: Decays widths and branching ratios of CP-odd A in the mssm model as a function of m^o, 
for the parameters M SU sy = 500 GeV, M 2 = 170 GeV, fj, = -A t = 1 TeV, tan (3 = 2.7. 
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FIG. 8: Widths decays and branching ratios of Higgs boson CP-odd ^4° in the MSSM model as 
function m A o, for the parameters Msusy = 500 GeV, M 2 = 170 GeV, \i = —A t = 1 TeV, 
tan = 20. 



IV. CONCLUSIONS 

We have calculated the one-loop contribution to the A°WW and A°ZZ couplings (absent 
at tree level) both in the context of Minimal Supersymmetric Standard Model MSSM and 
Two Higgs Doublet Model 2hdm. In the 2hdm, for small tan/3, near ti threshold, the 
branching ratio of A — > W + W~ and, A ZZ are enhanced to the level of 10 _1 and 10~ 2 
respectively. We stress also that in 2hdm type II and near ti threshold region, the branching 
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ratio of A — > gg could dominate aver A — > bb mode and hence A becomes fermiophobic. 
In the MSSM, the branching ratios of A -> W + W~ and A -> ZZ are of the order 10 
and 10 -4 respectively for low tan/3 and decrease for large tan/3. Those branching ratios, at 
this level, might provide an opportunity to search for a CP-odd Higgs boson at the lhc. 
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APPENDIX A: ONE-LOOP AMPLITUDES FOR A -» VV 

Let us briefly recall the definitions of scalar and tensor integrals we use. We use the 
convention of Loop Tools j^, 0|. The inverse of the propagator are denoted by, 

d = q 2 -ml , di = (q + pi) 2 - m\ (Al) 

where the pi are the momenta of the external particles (always incoming). 



in* J do"i 



Using Lorentz covariance, one gets for the vector integral 

B, = Pi M £i (A3) 

with the scalar function Bi(pl,mQ,m1). 
For the three-point function we have, 

c , = f (A4 ) 

«7r 2 J d did 2 

where p\ 2 = (pi +p 2 ) 2 - Lorentz covariance yields the decomposition 

C M = v\vP\ + P2^C 2 (A5) 
with the scalar functions Ci{j^,p\ 2 ,p\,ml,m\,m\). 
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APPENDIX B: LAGRANGIAN AND COUPLINGS 

In this Appendix, we use the fermion- vector boson coupling constants as defined in terms 
of the neutral-current and charged- current vertices. 

£ = ^w; |^7 m (^Pl + o$p R )x- + ^(Tl^Pl^ + vg KM d Lj yp L u Li ) + h.c} 

+ ^M^V(^ Pl + + ^(^JPl + Of Pa)xi- 

+ Yl /((- J / + Qf sin2 e w)P L + Q f sin 2 9 w P^j\ (Bl) 

f=l,u,d 

where the real 4 x 2-matrices Ofj R , symmetric and real 2 x 2-matrices 0'fj R , and symmetric 
and real 4 x 4-matrices 0"j R have the form: 

0% = N i2 U n - ^N i3 U j2 (B2) 

O r = N i2 V n + -j=NuV j2 (B3) 

0% = ~U i2 U j2 - UaUjx + sin 2 OwSti (B4) 

O y = -^-2^2 " VS-iKi + sin 2 0^- (B5) 



°"« = ^(NhN^-NoNx) (B6) 



2 



0"* = -(^t) (B7) 



The Lagrangian Higgs CP-odd neutralinos and charginos 

Ki,j=l i,j=l 

(B8) 

where 

Ofj NA = l - { ( sin /5iV 43 - cos {3N^ ( sin 6 w N n - cos W jV j2 ) + (i <- j) } (B9) 
5™ = (O™) (BIO) 
Og CA = -j= ( sin + cos /3C/ a ^ 2 ) (Bl 1) 



v/2' 

5goA = (B12) 
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APPENDIX C: ONE-LOOP AMPLITUDE 

In MSSM and 2HDM-II, the amplitude of the sum of diagrams Fig. ^a and ^b, is given 

by: 

Aw = A 2 V T + A^ Y (CI) 

• A — > ZZ 

^zz M = *4 E m W»f f { ({gSff) 2 + (9 L ZIf ) 2 ) C{ + (g* ff - g L zff ) *C( } (C2) 

Afr = E ^i^^ A {^^+(^+^j^+(^-^)<^3) 

- ^ E K ± (<<^ 



W i,j,fc=l 



-™ in' L n' L n CCA n' R n lR cCCA^^ikj 
- m^±\y jk u ik u^ - u jk u ik £> i:j jo 2 

Where fiU°/7 = tan /3 (cot f3), for up fermion (down fermion). The arguments of Cq X are 
(m%,m 2 z ,m 2 A0 ,mj,m 2 f ,m 2 f ) and C^( >2 are (m 2 z ,m 2 z ,m 2 A0 ,m%.,my k ,m%), where ^ a = x° or 



A -> W+W" 



= i E ^^//(CP + CT 7 ) (C4) 



yi 2susy _ .2mw" ^ n NNA$ ( n R n R , n L n L 

Cw ,jr, i v 



m^o(C + Ci) + (m^o - m^o)C 2 



(C5) 
(C6) 



fc=l,4 



where the Cq{ ^ have the same arguments Cq{ ^(m^, m^, m^o, m^, m^,, m^), all the Co,!^ 
have also the arguments Co^im^, m^r, m A o,m^,,m^,m^) where ^ = xf or x°. 
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• A -> Z 7 

^l| H 7 DM = -z2 tan % N c Qf ™f9A°ff (dzff + 9zjf) C ° ( C7 ) 

2 

^ U 7 SY = -^2m z sin^^(0'S + ^)(4 CA -Og CA ){m-± , C + (m-± - m-±) C 2 ) 
where Cq^ have the same arguments 6*0,2(0, m z , m\, m 2 .±, m1±, m 2 .±). 

Xi X{ Xj 

• A — > 77 

^2hdm = _ 24 £ iV c Q2 m 2 gin 2 e w g AOff C (C8) 

i=l,2 

The C have the argument C (0, 0, m^, m? ± , m?±, m?±). 
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